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Abstract We discuss the current bounds on the Equivalence Principle, in 
particular from structure formation and, reexamine in this context, the recent 
claim on the evidence of the interaction between dark matter and dark cn- 
. ergy in the Abell Cluster A586 and the ensued violation of the Equivalence 

Principle. 
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1 Introduction 



The Equivalence Principle (EP) is a basic physical ingredient of general rel- 
(f) ' ativity and shown to hold to great accuracy in composition dependent free 

fall experiments involving known forms of matter. However, recent discovery 
of empirical evidence of the existence of dark matter (DM) in the cluster of 
C*"** , galaxies 1.E0657 — 56, the so-called "bullet cluster" pQ, and in a ring-like DM 

structure in the core of the galaxy cluster C70024 + 17 [2], together with pro- 
posals on coupled dark energy (DE) have opened up the possibility of having 
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dark sector interactions which could alter the free fall of dark matter (DM) 
and DE as compared to baryons. 

In a recent paper [3] , we have argued that the study of the virial equilibrium 
of the cluster Abell A586 indicates DE-DM interaction and that this induces a 
detectable violation of the EP. This prompted a discussion on the connection 
with existing related work on the free fall in the presence of DM [H[6]as well 
as on the phenomenological bounds on the interaction of DE and DM arising 
from cosmology [7]. 

In this paper we examine the consistency of the results of Bertolami et 
al. [3] with the existing literature. Our discussion is organized as follows: in 
section[2]we briefly summarize the way EP violations can be detected. We then 
compare our previous results with homogeneous cosmology data in section [3] 
and express them in terms of free fall in structure formation in section |4] 
Section [5] contains our conclusions. 

2 An EP Violation Typology 

One can identify two qualitatively distinct forms of violation of the EP: i) A 
more classic one, where one observes a differential acceleration on the free fall 
of baryons in a gravitational field, the so-called fifth force type experiments 
and ii) A dark sector violation, where one compares the free fall of baryons 
with the one of DM or DE in a cosmological setting. 

2.1 Baryonic composition dependent violation of the EP 

In the fifth force type experiments one is concerned with the detection of 
the difference on the free fall acceleration of baryonic matter with different 
composition. A putative difference would be due to a composition dependent 
new interaction. Current bounds on such a differential acceleration imply that 
it should be smaller than 0.3 x 10~ 13 [5]. Alternatively, one could consider a 
differential gravitational acceleration for baryonic matter and for DM or DE 
(DM and DE if these two entities are interacting). The situation with DM- 
baryons composition dependent interaction was examined in Stubbs jT and the 
bound 3.7 x 10~ 3 for the differential acceleration of galactic DM on different 
baryonic test masses was obtained. The latest bound however is (—4 ± 7) x 
10~ 5 [S]. For sure, existing laboratory tests of the EP are irrelevant in what 
concerns its validity on the dark sector. For a general discussion of the EP 
validity for known particles and antiparticles, see, for instance, Bertolami et 
al. [9]. 

2.2 Dark sector violation of the EP 

The fact that recent observations suggest that about 96% of the content of the 
universe is composed by dark material of yet unknown nature opens up the 
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possibility of considering models where the free fall of dark sector particles is 
different from the one of the baryonic sector. Since the EP relies on clustering 
of mass, a violation of the EP would be detected by comparing the baryonic 
free fall with the clustering component of the dark sector: DM. That entails 
the two possibilities for such modification: either the self-gravity of DM gets 
modified in the particular fashion suggested in Kesden and Kamionkowski [5] , 
or this modification is induced by an additional interaction between clustered 
and unclustered components of the dark sector, DM and DE, respectively [3]. 
Contrary to the previous case, detection cannot rely on free fall experiments 
of different baryonic materials, but must involve baryon and DM falling onto 
DM 03]. 

Such configuration involves the formation of large scale structure, at ranges 
where DM plays a leading role. In Bertolami et al. [3], we have used the gravi- 
tational relaxation in a galaxy cluster, A586, and its peculiar type of statistical 
equipartition of energy to test for the presence of modifications with respect to 
the dust behaviour globally assigned to clustering matter. We have detected 
such modification and identified it with the DE-DM interaction, which itself 
logically entails a violation of the EP. We have proposed in Bertolami et al. [3] 
that such violation should be detectable in a generic bias present at the ho- 
mogeneous level. In section [4] we show that a particular model of DE-DM 
interaction yields a variation of the gravitational coupling equivalent to the 
effects of the ad hoc model considered in Kesden and Kamionkowski [5] . Such 
variation is used there to separate the free fall of baryons (stars) onto the 
DM filled galactic centre at a different rate than that of DM sub-haloes of the 
Milky Way. 



3 DE-DM Interaction and EP Violation 

Models of interactions in the dark sector have an impact on the EP since 
they modify the conservation of DM, at least at the homogeneous level. Such 
modification entails detectability on the evolution of the homogeneous geom- 
etry of the universe, as shown for the phenomenologically consistent unified 
model of DE and DM [10], the generalized Chaplygin gas (GCG) model [11] 
or for the more generic interacting model [Till 2] , but also on the quality of the 
self-gravitating equilibrium reached by large amounts of DM in clusters [3] . In 
Guo et al. [7], the consistency of two DE-DM interaction models, a constant 
coupling model and a varying coupling model [13], with respect to the con- 
junction of data from Cosmic Microwave Background (CMB) radiation shift 
parameter [14] , SNLS observations of type la Supernovae (SNIa) [15] and the 
Baryon Acoustic Oscillation (BAO) peaks in the SDSS [16]. The varying cou- 
pling model was also considered, together with the GCG model, in the study 
of the A586 cluster [3], in which detection of the DE-DM interaction was in- 
ferred from departure of the virial equilibrium. In what concerns the varying 
coupling model, observational evidence favors the detection of interaction in 
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Fig. 1 Superimposition of the probability contours for the interacting DE-DM model in the 
(<^x,£) plane (denoted as (lode ,v) m Bertolami et al. [3]), marginalized over QdEq in Guo 
et al. study of CMB, SNIa and BAO (2.66< £ <4.05 at 95% C.L.) with the extended 
results of Bertolami et al. [3] based on the study of A586 cluster. The £ = — 3lux dashed 
line corresponds to uncoupled models. The result of Bertolami et al. [3] is set at u)x = — 1 
and error bars are shown. The thick lines corresponds to their extension to varying u>x ■ 



the region where both approaches are compatible, namely for ojx < — 1 (see 
Fig. EE). 

The absence of detection in the homogeneous study [7] at 1 a level can 
be explained by the redshifts of CMB, SNIa and BAO, which are much larger 
than the one of A586, while in the regime of dominance of DE, the latter object 
is more affected. Should we have access to more systems like A586, an error 
analysis would be possible, adding the virial equilibrium as a Gaussian prior to 
obtain a joint likelihood analysis with the homogeneous constraints. With only 
one system to analyze, the signature shown on Fig. [1] merely superimposes the 
results from Guo et al. [7] with our study [3] with the respective errors. The 
overlap of contours from both studies favors interacting phantom-type models 
(ujx < —1) with energy transfer from DM to DE, as found in Bcnto et a/.and 
Majerotto et al. [TT1IT5]. 
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4 EP Violation, Free Fall and Structure Formation 

In their simulation of a particular modification of gravity for DM, Kesden 
and Kamionkowski [5] conclude that the remnants of stars (baryonic matter) 
in the heading tidal streams of the Sagittarius Milky Way satellite limits the 
variations induced by the resulting variation of gravitational coupling, G, up to 
10%. Considering a model for the DE-DM interaction together with constant 
mass DM particles, the effect on free fall gravitational force is also to modify 
the coupling with time: the DM particle number density, n, is diluted due to 
expansion and due to interaction with DE. This interaction can be modeled 
as an extra acceleration a/„t. felt by DM particles on top of a grav ., the gravity 
of the cluster. Since this extra acceleration is only felt by DM particles, the 
violation of the EP is explicit in this model. Furthermore, using that the 
relative change of total acceleration caused by this interaction is due to the 
relative change in particle number density with respect to expansion over a 
given cosmological time interval At: 

O-lnt. _ (P'DM — flDM Wexpansion) At 



oc . (1) 

a grav. n DM 

Since this is a cosmological effect, it is natural to consider At oc 1/H, the 
Hubble time. The constant DM mass hypothesis leads to pdm = riDMiriDM, 
which from the Bianchi identity, Eq. (1) from Bertolami et al. S\, turns into 

nDM + SHtidm = QHtidm , (2) 

where 

^ = (v + 3ujde)^de ^ 

Now, absence of interaction (£ = 0), Eq.0 would yield n DM \\ expanslon — 
SHnrjM- Thus, introducing a phcnomcnological constant 5 that also absorbs 
the Hubble time fraction, and using Eq. @, one obtains 

dint. ? n DM — n DAl\\ expansion rA /.s 

= ° T7 = 6 C ■ (4) 

Igrav. TlDMn 

Since this ratio depends only on the scale factor, a, the relative change in 
acceleration is conveyed directly to the effective DM gravitational coupling, 
Gdm ■ 

Gdm = G(l + <5C) , (5) 

which allows for a straightforward comparison with the simulation of Kesden 
and Kamionkowski 5J. As this change is due to the interaction with DE, 
which is significant only quite recently in the history of the universe, Gdm 
decreases as the strength of interaction becomes more important. Notice that 
the interaction vanishes in the past and then Gdm(z) = G. Given that the 
study of Kesden and Kamionkowski [S] sets from the observations a limit of 
about £ = 10% as the tolerance for such a change, we choose it to be the value 
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for Gdm at present, therefore in the present model 8 = e/|C| = 0.163, assuming 
that ludEo = ~ 1) ^de„ — 0.72, QdMq — 0.24, from which follows thau 
r) = 3.82^Q 47. Notice that the negative sign, given by our observation of A586's 
virial ratio, arises from the fact that the detected energy transfer is from DM to 
DE, i.e. one has less DM and hence effectively a smaller effective gravitational 
coupling. For sure, one does not expect a negative overall coupling as the 
discussed effect is presumably small and ours is just a phenomenological model 
of the complex process of interaction between concentrations of DM, which 
are at the same time turning into DE. It is interesting to remark that as the 
flow of energy density is from DM to DE, then udm is further decreased with 
respect to expansion, which translates into a decrease on the coupling of DM at 
present. Thus, one expects streams of stars to head the tidal arms of absorbed 
satellites in this model. 

Notice also that this calculation explicitly aims to translate the interaction 
with DE into a phenomenological variation of gravity in clusters. Of course, 
the value of e is chosen conservatively so that the modification is compatible 
with the findings of Kesden and Kamionkowski [5], However, it is fair to say 
that only a simulation with varying Gdm would be fully conclusive. Actually, 
it is quite conceivable that observations are consistent with even higher values 
of e. 

The behaviour of Eq. ([5]) is shown in Fig. [2] 

In the following, we discuss some of the salient features of the proposed 
model. For the Milky Way satellite, the observed stream evolution runs through 
at least 4 orbits that can be evaluated at 0.85 Gyr each. This corresponds to 
redshifts between about z = 0.36, where it started with Gdm/G = 0.94, and 
z — 0, with a midpoint in evolution time at z — 0.15 where the ratio is 0.91. 
For the A586, z = 0.17, and one obtains changes of Gdm of about 8% for the 
values of Bertolami et al. [3] (see above), as Eq. §5§ yields Gdm/G = 0.92. 
Notice that the observations of tidal streams from satellites of galaxies thus 
far should reveal the kind of trailing studied in Kesden and Kamionkowski [5] . 
[Furthermore, the negative contribution to the effective gravitational coupling 
can be attributed to the phantom-like nature of the interacting DE model.] 

It is relevant to point out that, in the context of the bullet cluster [z = 
0.296), an analysis of the velocity of the "bullet" sub-cluster indicates that a 
long range force with strength Gdm/G = 1.4 — 2.2 is needed to account for 
the rather high observed velocities |BJ, even though a more elaborate analysis 
[T9] suggests that more detailed hydrodynamical models are required to reach 
a definite conclusion. Our fit (which yields Gdm/G = 0.93, for z = 0.296 
with the values from Bertolami et al. [3]) does not account for such a large 
deviation from the Newtonian regime. And clearly, a 5 > is not even qualita- 
tively consistent with it, however, we stress that the parameter 6 reflects our 

1 We point out that errors in Bertolami et al. [3] were underestimated. The correct results 
are as follows: the virial ratio = —0.76 ± 0.14 7^ —0.50, the ui^e = — 1 interacting 
pseudo-quintessence scaling r\ = 3.82^JJ'^g ^ — 3ui£,E and the generalized Chaplygin gas 
power a = 0.27^q'1I 7^ . Conclusions remained unchanged. 



7 




Fig. 2 Evolution with redshift of the ratio of the gravitational coupling for DM and baryons 
falling on a DM halo, using the varying coupling model discussed in Bertolami et al. [3], as 
compared with the simulation of Kesden and Kamionkowski [5] . 



ignorance on the correct model of the interaction process and its details. In 
this context, it is interesting to point out that if instead, 5 were negative, then 
our approach could, at least in principle, be consistent with values required to 
accommodate the bullet cluster dynamics. 

Finally, it is worth remarking that our results are consistent with limits on 
the variation of the gravitational coupling with the energy scale |22] . 



5 Discussion and Conclusions 

In this paper we have examined the compatibility of the existing studies on 
the interaction of dark matter and dark energy and the results of Bertolami 
et al. [3] concerning the Abell A586 cluster. We have shown that there is no 
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contradiction with the available cosmological data, at least in what concerns 
the varying coupling model of interaction between DE and DM. Furthermore, 
we find that if one interprets the evidence of EP violation from the Abell A586 
as a change in the gravitational coupling, then our results match those obtained 
in Kesden and Kamionkowski [5] through a simulation of the tidal stream of 
stars in the Sagittarius dwarf galaxy due to DM. For sure, a more accurate 
comparison would involve considering a simulation carried out in Kesden and 
Kamionkowski 5J with the gravitational coupling varying according to Eq. 
However, this is clearly beyond the scope of this work. 
It should be mentioned that an attractive new force or a stronger gravi- 
tational coupling for DM might be an interesting solution to many small dis- 
crepancies one encounters when attempting to fit data with the ylCDM model. 
These discrepancies include, the greater than expected number of superclus- 
ters in the SDSS [20] , the extreme low density of matter within the voids [21] , 
the number of satellite galaxies which is an order of magnitude smaller than 
predicted |23 a among others. For sure, one can envisage other fixes for these 
problems, such as, for instance, self-interacting dark matter (see e.g. Bento 
et al. [24] and references therein). Whatever solution, one concludes that ev- 
idence for a complex and interacting dark sector is mounting. It is then just 
logical to question whether a basic physical ingredient such as the EP may 
hold for this sector as well. 
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